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Abstract 

The H*Htt form factor for H = B and D mesons is evaluated in a QCD sum 
\Q \ rule calculation. We study the Borel sum rule for the three point function 

of two pseudoscalar and one vector meson currents up to order four in the 
■ operator product expansion. The double Borel transform is performed with 

respect to the heavy meson momenta. We discuss the momentum dependence 
. of the form factors and two different approaches to extract the H*Hir cou- 

pling constant. 

Ph! PACS numbers 14.40.Lb, 14.40.Nd, 12.38.Lg, 11.55.Hx 
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.£h ! The coupling of the pion to the heavy mesons (gB*B-w and gD*D-w) is related to the form 

factor at zero pionic momentum and its precise value has been often needed in phenomenol- 
ogy. In particular, the Qo'D-k coupling is needed in the context of quark gluon plasma 
(QGP) physics. Suppression of charmonium production in heavy ion collisions is one of 
the signatures of QGP formation [0]. Therefore a precise evaluation of the background, 
i.e., conventional J/ip absorption by co- moving pions and p mesons 0, is of fundamental 
importance. Since pions are so abundant in a dense nuclear environment, the reactions 
ir + J/ip^D + D* (and consequently the coupling gD*Dx) are of special relevance [[|. 

In the case of gD*Diz, the D* + — > D°tt + decay is observed experimentally. However, 
present data provide only an upper bound: gD*Dn < 21 |§]. For gB*Bn, there cannot be 
a direct experimental indication because there is no phase space for the B* — ► Bn decay. 
Recently, a direct preliminary determination of gB*B-K on the lattice has been attempted ||. 

The D*Dn and B*Bn couplings have been studied by several authors using different 
approaches of the QCD sum rules (QCDSR): two point function combined with soft pion 
techniques 0,0], light cone sum rules light cone sum rules including perturbative correc- 
tions ||10|| , sum rules in a external field pT| , double momentum sum rules []T2| | . Unfortunately, 



the numerical results from these calculations may differ by almost a factor two. 
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In this work we use the three-point function approach to evaluate the D*Dir and B*Bir 
form factors and coupling constants. The advantage of using the three-point function ap- 
proach with a double Borel transformation compared with the two-point function with a sin- 
gle Borel transformation is the elimination of the terms associated with the pole-continuum 



transitions P|,|T3"lj. 



The three-point function associated with a H*Hn vertex, where H and H* are respec- 
tively the lowest pseudoscalar and vector heavy mesons, is given by 

r>,p') = / d 4 x A(0|T{j(x)j 5 (l/)jJ(0)}|0)e^e-*'-^^ , (1) 

where j = iQ^u, j$ = iu'j^d and ft = dj^Q are the interpolating fields for H, ir~ and H* 
respectively with u, d and Q being the up, down, and heavy quark fields. 

The phenomeno logical side of the vertex function, I^Qo, p'), is obtained by the consider- 
ation of H and H* state contribution to the matrix element in Eq. ([]]): 

rf en \p, P ') = 2 1 2 /2 1 2 (o\j\Htf)) x 

p p z — m H * p' — m H 

(H(p)\j 5 \H*(p, e)){H*(p, 6)1^10) + higher resonances . (2) 

The matrix element of the pseudoscalar element, js, defines the vertex form factor 
9h*hM 2 )- 



(H(p')\j 5 \H*(p,e)) = £g* 9H ' 2 HM l\ ^ , (3) 

m u + m d q z - m z 



where q = p' —p, f„ is the pion decay constant and e y is the polarization of the vector meson. 
The vacuum to meson transition amplitudes appearing in Eq. (0) are given in terms of the 
corresponding meson decay constants fn and fn* by 

(0\j\H(p')) = ^ , (4) 

m Q 

and 

(H*(p,e)\jl\0) = m H *f H *e;. (5) 
Therefore, using Eqs. (|3|), ([|) and (||) in Eq. (g) we get 

rr eri W) = c hh j h ; h ^ 2 1 2 , 2 1 2 x 

p q z — m z p z — m H * p' — m H 

( , m 2 H * + m 2 H — q 2 \ 

—p., H 5 Pa + higher resonances , (6) 

y ^ 2m H * I 



where 



c ^ = m 2 H m H *mlf H f H *f n 
HH (m u + m d )m Q 
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The contribution of higher resonances and continuum in Eq. (§) will be taken into account 



as usual in the standard form of ref. 14 



The QCD side, or theoretical side, of the vertex function is evaluated by performing 
Wilson's operator product expansion (OPE) of the operator in Eq. (|]). Writing in terms 
of the invariant amplitudes: 



^(p, p) = Ti(p 2 , p' 2 , q 2 )p^ + T 2 (p 2 , p'\ q 2 )p[ 



,2 



we can write a double dispersion relation for each one of the invariant amplitudes Tj (i 
over the virtualities p 2 and p' 2 holding Q 2 = —q 2 fixed: 



Up 2 ,p' 2 ,Q 2 



4vr 2 



ds 



p 2 ){u — p' 2 ) 



(8) 
1,2), 

(9) 



where pi(s, u, Q 2 ) equals the double discontinuity of the amplitude Ti(p 2 ,p' 2 , Q 2 ) on the cuts 
ttiq < s < 00, nig < u < 00, which can be evaluated using Cutkosky's rules |I4| , |I5| . 

Finally we perform a double Borel transformation [|14[] in both variables P 2 = —p 2 and 
P' 2 = —p' 2 and equate the two representations described above. We get one sum rule for 
each invariant function. In the p^ structure: 



- C 



HH* 



m H* +m 2 H + Q 2 g H *HM 2 ) 

Q2 _|_ m 2 



™ 2 H JM\-ml/M'- 



1 



2m 2 H , 



ds 



du 



4tt 2 

Pl (s,u,Q 2 )e- s / M2 e- u / M ' 



(10) 



and in the p'^ structure: 
9H*H7r(q 2 ) 



C 



HH* ' 



Q2 _)_ m 2 



m%»/M 2 -m 2 H /M' 2 



4vr 2 



Uo 



ds / du 



p 2 (s,u,Q 2 )e- s / M2 e- u / M ' 



(11) 



where So and Uq are the continuum thresholds for the H* and H mesons respectively, which 
are, in general, taken from the mass sum rules. The two Borel masses M 2 and M' 2 are, 
in principle, independent and they should vary in the vicinity of the corresponding meson 
masses: m 2 H * and m 2 H respectively. Since for heavy mesons m# and m#* are very close, 
many authors use M 2 = M' 2 To allow for different values of M 2 and M' 2 we take 



them proportional to the respective meson masses, which leads us to study the sum rule as 
a function of M 2 at a fixed ratio 



M 2 



m 



H* 



ni 



(12) 



H 



We will consider diagrams up to dimension four which include the perturbative diagram 
and the gluon condensate. The quark condensate term does not contribute since it depends 
only on one external momentum and, therefore, it is eliminated by the double Borel transfor- 
mation. Higher dimension condensates are strongly suppressed in the case of heavy quarks 
p|-p|, pd] , p!2| . The double discontinuity of the perturbative contribution reads: 
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(13) 



{vert), n2 , SQ 2 [m 2 Q (s + u + Q 2 ) -2su] 
* (g -"' g)= [(s + u + Q 2 ) 2 -Asu}V 2 ' (14) 

and the integration limit condition is 

(s-m 2 Q )(u-m 2 Q )>Q 2 m 2 Q . (15) 

In this paper we focus on the structure which we found to be the more stable one. 
For consistency we use in our analysis the QCDSR expressions for the decay constants up 
to dimension four in lowest order of a s : 



= f: - £msw- . tie) 

o7T iiljj JtDp U ''^H 



ft. = f S ° ds {S ~ m2Q? (2 + ^) e^-V^ 9 - ^-(qq)e^*- m , (17) 

8n 2 m z H * Jm? Q s \ s J m z H * 

where we have omitted the numerically insignificant contribution of the gluon condensate. 

The parameter values used in all calculations are m u + = 14 MeV, m c = 1.5 GeV, 
m b = 4.7 GeV, m„ = 140 MeV, m D = 1.87 GeV, m D * = 2.01 GeV, m B = 5.28 GeV, 
m B , = 5.33 GeV, / w = 131.5 MeV, (qq) = -(0.23) 3 GeV 3 , (g 2 G 2 ) = 0.5 GeV 4 . We 
parametrize the continuum thresholds as 

so = {m H * + A s ) 2 , (18) 

and 

u = (rnji + A u ) 2 . (19) 

The values of uq and so are, in general, extracted from the two-point function sum rules for 
f H and f H * in Eqs. (||) and (0). Using the Borel region 2 < M 2 < 5 GeV 2 (for the D* and 
D mesons) and 10 < M 2 < 25 GeV 2 (for the B*, and B mesons) we found a good stability 
for fjj and fn* with A s = A u ~ 0.5 GeV, in agreement with the results in ref. ||. We have 
checked that bigger values for A s ( u ), of order of 1 GeV, lead to unstable results for f H and 
fn*, in the case of the sum rules Eqs. fll6|) and (p7T|) . In our study we will allow for a small 
variation in A s and A u to test the sensitivity of our results to the continuum contribution. 

We first discuss the D*Dir form factor. In Fig. 1 we show the behavior of the perturbative 
and gluon condensate contributions to the form factor gn*Dn(Q 2 ) at Q 2 = 1 GeV as a 
function of the Borel mass M 2 using A s and A n given in Eqs. ( [T8] ) and flT9"|) equal to 0.5 GeV. 
We can see that, in the case of the form factor, the gluon condensate is not negligible and 
it helps the stability of the curve, as a function of M 2 , providing a rather stable plateau 
for M 2 > 3 GeV 2 . The behavior of the curve for other Q 2 and continuum treshold values 
is similar. Fixing M 2 = 3.5 GeV 2 we show, in Fig. 2, the momentum dependence of the 
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form factor (dots). Since the present approach cannot be used at Q 2 = 0, to extract the 
9d*Dtt coupling from the form factor we need to extrapolate the curve to Q 2 = (in the 
approximation m 2 =0). In order to do this extrapolation we fit the QCD sum rule results 
(dots) with an analytical expression. We tried to fit our results with a monopole form, since 
this is very often used for form factors, but the fit is very poor. We obtained good fits using 
both the gaussian form 



9H*Htt(Q 2 ) — QH'H-kZ 



-(Q 2 + m 2)Vr 4 



and a curve of the form 



gH*Hir(Q 2 ) — 9h*h-k 



1 + {a/A) 4 



1 + (a/A)W+ m ') 2 /A 4 



(20) 



(21) 



In Fig. 2 we show that the Q 2 dependence of the form factor, represented by the dots, can be 
well reproduced by the parametrization in Eqs. ( ^0|) (dashed line) and (^) (solid line). The 
value of the parameters in Eqs. (^0|) and (^1]) are given in Table I for two different values of 
the continuum threshold. 



A s = A u ( 


GeV) 


9D*Dir 


A ( GeV) 


a ( GeV) 


T ( GeV) 


0.5 




5.3 


1.66 


1.90 




0.6 




6.0 


1.89 


3.05 




0.5 




5.7 






1.74 


0.6 




6.1 






1.92 



TABLE I: Values of the parameters in Eqs. (^) and (^) which reproduce the QCDSR results 
for 5d*Dtt(Q 2 )) for two different values of the continuum thresholds in Eqs. ( |l8l) and dill). 



In view of the uncertainties involved, the results obtained with the two parametrizations are 
consistent with each other, the systematic error being of the order of 10%. 

In refs. P,l~6| it was found that the form factor in the semileptonic decay H — > ttW, which 
is also normalized by the H*Hir coupling constant, can be well approximated by a monopole 
form factor. In the case of the H — > nlv form factor, a vector dominance approximation 
gives a phenomenological explanation for a pole fit at q 2 = m^*, which is not the case of 
the form factor studied here. It is important to notice that here the dispersion relation is 
written in terms of the two heavy meson momenta, while in the case of semileptonic decay 
the dispersion relation is a function of the H and 7r momenta. Therefore, our form factor is 
a function of the pion momentum, exhibiting a peak at the pion pole Q 2 = 0. 

To test if our fit gives a good extrapolation to Q 2 = we can write a sum rule, based 
on the three-point function Eq. (|1|), but valid only at Q 2 = 0, as suggested in |T7| for the 
pion-nucleon coupling constant. This method was also applied to the nucleon-hyperon-kaon 
coupling constant [|T^,|T^] and to the nucleon— A c — D coupling constant ||20[| . It consists in 
neglecting the pion mass in the denominator of Eq. (§) and working at Q 2 = 0, making a 
single Borel transformation to both P 2 = P' 2 — > M 2 . 

As discussed in the introduction, the problem of doing a single Borel transformation is 
the fact that the single pole contribution, associated with the N — > N* transition, is not 
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suppressed |||]|L3). In ref. |T3[ it was explicitly shown that the pole-continuum transition 
has a different behavior as a function of the Borel mass as compared with the double pole 
contribution and continuum contribution: it grows with M 2 as compared with the double 
pole contribution. Therefore, the single pole contribution can be taken into account through 
the introduction of a parameter A, in the phenomenological side of the sum rule P, |i~3| , |T9" 



Thus, neglecting m 2 in the denominator of Eq. @ and doing a single Borel transform in 
P 2 = P' 2 , we get for the structure p^ 

fr n) (M 2 , Q 2 ) = _ C H , H m 2 H + m 2 H , + Q 2 , m y M , _ ^jm^ ( + ^ 

2m z H Q 2 m z H , — m A H v ' 

(22) 

where Ch*h in given in Eq. (^) with f H and fn* given by Eqs. (|1^) and (|17D. 

On the OPE side only terms proportional to 1/Q 2 will contribute to the sum rule. 
Therefore, up to dimension four the only diagram that contributes is the quark condensate 
given by 

f<^>(M 2 , Q 2 ) = e -™%/M> (23) 

Equating Eqs. (|22| ) and (|23| ) and taking Q 2 = we obtain the sum rule for gH*Hn + AM 2 , 
where A denotes the contribution from the unknown single poles terms. It is interesting to 
point out that in the limit m 2 H + m 2 H * = 2m|p , the sum rule obtained in the p'^ structure 
coincides with the sum rule in the p^ structure. In Fig. 3 we show, for A s = A u = 0.5 GeV, 
the QCDSR results for go'D-w + AM 2 as a function of M 2 (dots) from where we see that, in 
the Borel region 2 < M 2 < 5 GeV 2 , they follow a straight line. The value of the coupling 
constant is obtained by the extrapolation of the line to M 2 = | I3| . Fitting the QCDSR 
results to a straight line we get 

9D*Dn ^ 5.4 , (24) 

in excellent agreement with the values obtained with the extrapolation of the form factor to 
Q 2 = 0, given in Table I. 

It is reassuring that both methods, with completely different OPE sides and Borel trans- 
formation approaches, give the same value for the coupling constant. 

In the case of B*E>7r vertex, we show in Fig. 4, for A s = A u = 0.5 GeV, the Q 2 = sum 
rule results for gB*Bn + AM 2 (dots) as a function of M 2 . It also follows a straight line in 
the Borel region 10 < M 2 < 25 GeV 2 , and the extrapolation to M 2 = gives 

9B* Bit ^ 10.6 . (25) 

In Fig. 6 we show the QCDSR result for the perturbative and gluon condensate con- 
tributions to the form factor gB*Bw(Q 2 ) at Q 2 = 2 GeV as a function of M 2 using 
A s = A u = 0.5 GeV. In this case the gluon condensate is very small but it still goes in the 
right direction of providing a stable plateau for M 2 > 15 GeV 2 . Fixing M 2 = 17 GeV 2 we 
show, in Fig. 6, the Q 2 behavior of the form factor (dots). The dots can still be well fitted 
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by Eq. ([H]) (solid line). However, the fit with Eq. ( p0| ) is not so good, as can be seen by the 
dashed line in Fig. 6. In Table II we give the value of the parameters in Eqs. (|2(|) and ( |2l"l) 
that reproduce our results for two different choices of the continuum thresholds. 

In this case the agreement of the two different approaches to extract the coupling constant 
is not so good, but the numbers are still compatible. One possible reason for that is the 
fact that for heavier quarks the perturbative contribution (or hard physics) becomes more 
important, as can be observed by the decrease of the importance of the gluon condensate in 
Fig. 5 as compared with Fig. 1. 



Since in the sum rule given by Eqs. (|22|) and (|23"D there is 
only soft physics information, we expect a s corrections to the sum rule to be more important 
in the case of gB*Biv{Q 2 ) than for 5 , d*Dtt(<5 2 )- 



A s = A u ( GeV) 


9B* Bit 


A ( GeV) 


a ( GeV) 


T ( GeV) 


0.5 


14.7 


1.62 


1.37 




0.6 


16.3 


1.81 


1.67 




0.5 


17.2 






1.79 


0.6 


18.4 






1.97 



TABLE II: Values of the parameters in Eqs. (20) and fl2T|) which reproduce the QCDSR results 
for gB*BTr(Q 2 ), for two different values of the continuum thresholds in Eqs. ( |T8| ) and (|l9|). 



Comparing Table I with Table II we see that the cut-offs are of the same order in the two 
vertices and are very hard. Concerning the parameter a, it is smaller in the case of the 
B*Bn vertex. This is because of the fact that the form factor gB*Biv{Q 2 ) has a flatter peak 
around Q 2 = than gD*Dir{Q 2 )- This can be interpreted as an indication that the spatial 
extension of the vertex is smaller for B*Bir than for D*Dn. This is also the reason why 
the gaussian fit is not so good in the case of the B*Bti vertex, and leads to bigger values 
for the coupling. It is interesting to notice that our results for the coupling constants are 



completely consistent with the QCDSR calculation of ref. |L2 . 

As a final exercise, we use our result for gs*B-n to extract the coupling constant g which 
controls the interaction of the pion with infinitely heavy fields in effective lagrangian ap- 
proaches [gT],^[. They are related by [0-|,|TJ,|rj,fn],£g] 

2mB (oa\ 
gB*B-K = —^9 ■ (26) 

The knowledge of g is of great phenomeno logical value, since its strenght is required in the 
analyzes of many electroweak processes [[Hp. Therefore, during the last years, a large number 
of theoretical papers has been devoted to the calculation of g. However, the variation of 
the value obtained for g, even within a single class of models, turns out to be quite large. 



For instance, using different quark models one obtains 1/3 < g < 1 [22,23] while QCDSR 
calculations points in the direction of small g, with a typical value in the range g ~ 0.13—0.35 

Mai. 



Using the values for gs*B-K given in Table II we get, at order a s = 

# = 0.13-0.23, (27) 
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therefore, we corroborate the overall conclusion drawn from different QCDSR calculations, 
that the coupling g is small. 

In conclusion, we extracted the H*Htt coupling constant using two different approaches of 
the QCDSR based on the three-point function. We have obtained for the coupling constants: 

<7 D * D7r = 5.7±0.4, (28) 



g B * Bn = 14.5 ± 3.9 , (29) 

where the errors reflect variations in the continuum thresholds, different parametrizations 
of the form factors and the use of two different sum rules. There are still sources of errors in 
the values of the condensates and in the choice of the Borel mass to extract the form factor, 
which were not considered here. Therefore, the errors quoted are probably underestimated. 

In Table III we present a compilation of the estimates of the coupling constants gD*Dir 
and gs*B-K from distinct QCDSR calculations. 



approach 


9d*d-k 


9B*Btt 


this work 

two-point function + soft pion techniques (2PFSP) 
2PFSP + perturbative corrections || 

light cone sum rules (LCSR) || 
LCSR + perturbative corrections [10] 
double momentum sum rule l2fl 


5.7 ±0.4 
9 ± 2 
7±2 
11 ± 2 
10.5 ±3 
6.3 ±1.9 


14.5 ±3.9 
20 ±4 
15 ±4 
28 ±6 
22 ±9 
14 ±4 



TABLE III: Summary of QCDSR estimates for gn^D-K and gB*B 



From this Table we see that our result is in a fair agreement with the calculations in 
refs. PJT^] , while LCSR calculations point to bigger values for the coupling constants. This 
discrepancy has still to be solved. 

The D*D7r coupling is directly related with the D* — > Dti decay width through 

2 I — *■ 1 3 

rm*- _> Wn-) = 9D * D ^y . (30) 

247rm|,« 

Using Eq. ( ^8l) we get 

T(D*~ -> Wn-) = 6.3 ± 0.9 keV , (31) 
which is much smaller then the current upper limit T(D*~ — > D°tt~) < 89 keV. 
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FIG. 1. M 2 dependence of the perturbative (long-dashed line) and gluon condensate (dashed 
line) contributions to the D*Dir form factor at Q 2 = 1 GeV 2 (solid line) for A s = A u = 0.5 GeV. 

7 i ' 1 1 1 1 1 1 1 1 1 1 1 
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1 2 3 4 5 6 
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FIG. 2. Momentum dependence of the D*Dtt form factor for A s = A u = 0.5 GeV (dots). The 
solid and dashed lines give the parametrization of the QCDSR results through Eqs. and (|20| ) 
respectively. 
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M 2 (GeV 2 ) 



FIG. 3. D*Dir coupling constant as a function of the squared Borel mass M 2 from the QCDSR 
valid at Q 2 = (dots). The straight line gives the extrapolation to M 2 = 0. 




M 2 (GeV 2 ) 



FIG. 4. B*Bir coupling constant as a function of the squared Borel mass M 2 from the QCDSR 
valid at Q 2 = (dots). The straight line gives the extrapolation to M 2 = 0. 
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FIG. 5. M 2 dependence of the perturbative (long-dashed line) and gluon condensate (dashed 
line) contributions to the B*Btt form factor at Q 2 = 2 GeV 2 (solid line) for A s = A u = 0.5 GeV. 
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1 2 3 4 5 6 

Q 2 (GeV 2 ) 

FIG. 6. Momentum dependence of the B*Btt form factor for A s = A u = 0.5 GeV (dots). The 
solid and dashed lines give the parametrization of the QCDSR results through Eqs. ( ^Tj) and (|20|) 
respectively. 
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